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1 Introduction 


If CP conservation holds in the system, the CP-even and CP-odd eigenstates are mass eigen¬ 
states with masses M+, M_ and widths r_|_, r_. The mixing parameters x = 2(M_|_—M_)/(r+-|-r_) 
and y = (r+ — r_)/(r+ -|- r_) measure the difference of these masses and lifetimes, respectively. 
Since in the Standard Model D-D mixing is doubly Cabibbo-suppressed and vanishes in the SU (3) 
flavor limit both parameters are expected to be small (10“^) |]|]. If D - D mixing were found to 
be large, this might be a sign either for large SU (3) flavor breaking (y > x) or for new physics 
(x > y) H- 

The mixing parameter y can be determined from the lifetime of D^ meson^ that decay into final 
states of specific CP symmetry [^. A final state that is an equal mixture of CP-even and CP-odd 
is produced by the Cabibbo-favored decay D^ —> K~'k~^ . If y is small, the lifetime distribution 
of D^ mesons decaying into this final state can be approximated as an exponential with lifetime 
tk-k = 1/r where T = (r+ -|- r_)/2. A CP-even final state is produced by the singly Cabibbo- 
suppressed decay D^ —> K~K^. The decay time distribution of D^ mesons that decay into K~K~^ 
is exponential with a lifetime tkk = l/r+- This lifetime can be compared to tk-k to obtain y: 


y = 


Tkk 


- 1 . 


( 1 ) 


The mixing parameters x^ and y can be determined simultaneously from the time evolution of 
the wrong-sign (WS) decay D^ —> K~^n~. The WS decay rate has contributions from the doubly 
Cabibbo-suppressed decay, described by a pure exponential, and from D-D mixing, described 
by the same exponential modified by a coefficient quadratic in the decay time t. A third term 
arises from the interference of the two, where the exponential has a coefficient linear in t. The time 
evolution of the WS decay rate is thus described by [P: 

r(t) oc exp(— t) [ R -|- \/p y' t 1/4 ( x'‘^ +y'‘^ ) ]• (2) 


Here, t is given in units of the lifetime of the D^. The parameters x' and y' are related to the 
mixing parameters x, y by a rotation: 


x' = X cos 6 + y sin 6, y' = y cos 6 — x sin 5. 


(3) 


The phase of the rotation is the strong phase between the doubly Cabibbo-suppressed contribution 
and the one from mixing. 

In the following, two results are presented. Firstly, a preliminary measurement is described 
of y from the difference of the D^ lifetime determined in the and K~K~^ decay channels. 

Secondly, a preliminary measurement of the WS decay rate of the D^ is discussed, extracted from 
a full fit to the time evolution of the combined right-sign and wrong-sign D^ decay rates. 

The measurements are based on data collected with the BABAR detector at the PEP-II asym¬ 
metric e'’'e“ collider. Data taken on and off the T(4S') resonance are used. Their center-of-mass is 
boosted along the beam axis with a nominal Lorentz boost of (3^ = 0.56. The size of the interaction 
point (IP) transverse to the beam direction is typically 6 ym in the vertical direction and 120 ym 
in the horizontal direction. 

Candidates for D^ particles are identified through the decay D*~^ —> D^Tr~^. Charged particles 
are detected and their momenta are measured by a combination of a 40-layer drift chamber (DCH) 

^In the following, the charge conjugates of the meson and its decay modes are always implied as well. 


2 




and a five-layer, double-sided, silicon vertex tracker (SVT), both operating in a 1.5 T solenoidal 
magnetic field. A ring-imaging Cherenkov detector (DIRC) is used for charged particle identifica¬ 
tion. A detailed description of the BABAR detector is available in |^]. 

2 Measurement of y 

The result is obtained from a sample of 12.4 fb~^ of 2001 BABAR data that were reconstructed with 
the most advanced tracking alignment parameters and reconstruction algorithms. 

2.1 Event Selection 

The widths r_ and r_|_ are determined by fitting the decay time distributions of independent 
samples of K~tt^ and —> K~K~^ decays. The candidates for each sample are 

identified by means of the charged particles in their final state. The decay D*~^ —> D^Tr~^ and 
particle identification are used to suppress backgrounds. 

candidates are selected by searching for pairs of tracks of opposite charge and combined 
invariant mass near the mass mo- Each track is required to contain a minimum number of 
SVT and DCH hits in order to ensure reconstruction quality. The two daughter tracks of the 
candidate are fitted to a common vertex. The probability of this vertex fit has to be better than 
1 %. 

Each candidate among the daughter tracks is required to pass a likelihood-based particle 
identification algorithm. This algorithm is based on the measurement of the Cherenkov angle by 
the DIRC for momenta p > 0.6 CeV/c and on the energy loss (dE/dx) measured by SVT and DCH 
for momenta p < 0.6GeV/c. An average efficiency for identification of approximately 80% is 
reached for tracks within the DIRC acceptance while the tt^ misidentification probability amounts 
to about 2%. 

The decay D*~^ —> D^Tr~^ can be distinguished by a vr"*" of low momentum, commonly referred 
to as the slow pion (tTs). To increase acceptance, vr^ candidate tracks do not have to contain DCH 
hits. To improve momentum resolution, a vertex fit is used to constrain each vr^ candidate track to 
pass through the intersection of the trajectory and the IP. If the probability of this vertex 
fit is less than 1%, the D* candidate is discarded. 

The D* candidates peak at a value of 5m ~ 145.4MeV/c^, where 5m is the difference in the 
reconstructed D* and masses. Backgrounds are reduced by rejecting events with a value of 5m 
that deviates more than a given margin from the peak. The size of this margin corresponds to 
approximately three standard deviations and varies between I and 2.5MeV/c^, depending on the 
quality of the vr^ track. 

To remove background from B meson decays, each D* candidate has to have a momentum p* 
greater than 2.6GeV/c in the center-of-mass. This condition is also effective at removing combina¬ 
torial background that tends to accumulate at lower momenta. 

After all selection requirements, about 45,000 candidates in the decay channel K~7r~^ and 
about 4,000 in the K~K~^ decay channel are retained. The mass and 5m distribution of the selected 
events are shown in Fig. |l|. The relative size of the remaining background is about 2% and 5% for 
the K~tt^ and K~K~^ samples, respectively, when measured inside a ±20MeV/c^ mass window. 
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Figure 1: The reconstructed 6m and mass distributions after event selection for the and 

K~K~^ decay modes. The 6m plots include candidates both inside and outside the 6m selection 
requirement which fall within the m^, window indicated in the lower plots. 


2.2 Lifetime Determination 

The flight length and its measurement error are determined for each candidate by a global, 
three-dimensional, multiple-vertex fit that includes the daughter tracks, the vr^ track, and the 
IP envelope. This fit does not include explicit constraints on the D* or masses. The value listed 
by the Particle Data Group (PDG) for the mass {mjj = 1.8654 GeV/c^ |^]) and the momentum 
of the obtained with the vertex fit are used to calculate the boost of the and obtain the 
proper decay time. 

The lifetimes of the candidate samples are extracted by means of an unbinned maximum 
likelihood fit. The likelihood function is divided into two distinct decay time distribution func¬ 
tions, one for the signal and one for the background. The signal function is the convolution of an 
exponential and a resolution function. The resolution function consists of the sum of two Gaussian 
distributions with zero mean and with widths that are proportional to the measurement error (typ¬ 
ically 180 fs) of the decay time of each candidate. The parameters in the fit associated with 
the signal are the lifetime, the proportionality factors for the widths of the two Gaussians, and the 
fraction of signal that is assigned to the second Gaussian. 
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Figure 2: The fit to the reconstructed decay time distribution in the two decay modes and 
for all events including the mass sidebands. The white histogram represents the result of the 
unbinned maximum likelihood fit described in the text. The gray histogram is the portion assigned 
by the fit to background. 


Like the signal likelihood function, the background function consists of the convolution of a 
resolution function and a lifetime distribution. The background lifetime distribution is determined 
as the sum of an exponential distribution and a delta function at zero, the latter corresponding to 
those sources of background that originate at the IP. The resolution function consists of the sum 
of three Gaussian distributions. The first two are chosen to match the resolution function of the 
signal. The third is given a width independent of the decay time error and accounts for outliers 
produced by long-lived particles or reconstruction errors. The additional fit parameters associated 
with the background include the fraction assigned to zero lifetime sources, the background lifetime, 
the width of the third Gaussian, and the fraction of background assigned to the third Gaussian. 

The reconstructed mass of each candidate provides the likelihood of this candidate to be part 
of the signal. This likelihood is based on a separate fit of the reconstructed mass distribution. 
This fit includes a resolution function composed of a Gaussian with an asymmetric tail and a linear 
portion to describe the background. The slope of the background part is constrained with 
candidates in the 5m sideband (151 < 5m < 159MeV/c^). 

The results of the lifetime fits are shown in Fig. Typical values for the fitted parameters are 
a background lifetime similar to the lifetime and a width of the third Gaussian that is several 
times larger than the typical decay time error. The proportionality factors associated with the two 
Gaussians in the resolution function correspond to a root-mean-square of approximately 1.2. 
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To ensure that the analysis was performed in an objective manner, the lifetime and y values 
were blinded throughout the analysis process. This blinding was performed by adding to each of 
the tk-k and tkk ht results an offset chosen from a random Gaussian distribution of width 10 fs. 
The values of the two offsets and the positive (r* > 0) sides of the lifetime distribution of the data 
and of the fit result (Fig. were concealed. The value of y was unblinded only after the analysis 
method and systematic uncertainties were hnalized and the result was committed for public release. 

2.3 Systematic Errors and Result 

Many systematic uncertainties cancel because y is measured from the ratio of lifetimes. The 
uncertainties that do not cancel are associated mostly with backgrounds. These are tested by 
varying each event selection requirement within its uncertainty and recording the subsequent change 
Aj in the measured value of y. The quadrature difference (<5Aj)^ = Ictq — (t?| is used as an estimate 
of the statistical error (5Aj in A,, where gq {at) is the statistical error in y before (after) the ith 
systematic check. Only systematic checks with Aj > SAi are included in the sum Af — {6Ai)‘^. 
The square root of this sum (1.7%) provides an estimate of the systematic uncertainty originating 
from the event selection and backgrounds. 

Biases in tracking reconstruction are investigated by studying Monte Carlo samples, which, 
within statistics, show no reconstruction bias. In addition, a variety of vertexing techniques, includ¬ 
ing constraining the mass and using separate D* and vertex fits are employed as systematic 
checks. A systematic uncertainty of 0.4% is assigned from this source. 

Detector misalignment is another potential source of systematic uncertainty. Residual internal 
misalignment of the SVT, even as small as a few microns, can produce significant variations in the 
apparent lifetime. 

This source is studied with the help of e'*'e“ ^ 77 —> events in which the four charged 

tracks are known to originate from the IP. By selecting oppositely charged pairs of these tracks with 
opening angles similar to two-prong decays, it is possible to measure the apparent IP position 
as a function of trajectory and to calculate a correction to the lifetime obtained by the fit. 
For the data samples used in this analysis, a correction of -|-5 fs is found, with negligible statistical 
error and a systematic uncertainty of ±5 fs. This type of correction nearly cancels in the lifetime 
ratio and introduces little systematic uncertainty in y. 


Category 

Uncertainty (%) 

Event Selection and Background 

1.7 

Reconstruction and Vertexing 

0.4 

Alignment 

0.3 

Quadrature Sum 

1.7 


Table I: Individual contributions to the systematic uncertainty in y. 

The systematic uncertainties in y are summarized in Table ||. Adding the contributions of all 
systematic checks in quadrature, the preliminary result for y is obtained: 

2 /= (-1.0 ±2.2 ±1.7) % , (4) 

where the first error is statistical and the second systematic. This result is consistent with the 
Standard Model expectation of zero and with results from E79I |@], FOCUS Belle |p and 
CLEO 
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An important test of the soundness of the analysis is a lifetime that agrees with the PDG 
value of 412.6±2.8 fs |^]. A corrected value of tk-k = 412±2 fs is found. The systematic uncertainty 
is approximately 6 fs and is dominated by detector alignment effects. 

3 Measurement of the wrong-sign decay rate 

Wrong-sign (WS) decays of the type —> K~^'k~ are selected by requiring that the slow pion 

and the kaon have the same charge. The measurement is based on 23fb~^ of BaBar data taken 
in 2000. Vertex reconstruction, event selection and systematic error analysis are similar to the 
ones employed in the measurement of y, described above. After all event selection requirements a 
sample of about 200 WS candidate events is retained. 

A simultaneous unbinned maximum likelihood fit to the decay time distribution of right-sign 
(RS) and wrong-sign decay candidates is performed. The fit function consists of a lifetime 
distribution function convoluted with a resolution function. For RS events, a pure exponential is 
assumed as functional form of their lifetime distribution while for WS events the functional form 
follows Eq. 1^ The simultaneous fit of RS and WS candidates exploits the high statistics of the RS 
sample to obtain an accurate description of the resolution function and of the type and size of the 
background contributions. 

The 1 and 2 a contours of the statistical error of the fit result in the y'-x'^-plane (Fig. exhibit 
a strong correlation between the two parameters that extends into the unphysical region < 0. 
Note that a length scale is given for the axes but no axis labels and thus no central values for 
y'. The analysis of systematic errors shows that the dominating systematic error source in the 2000 
BABAR data is the internal alignment of the SVT. R will soon be possible to extract the central 



Figure 3: The 1 and 2 a contours of the statistical error in the x'^-y^-plane {y' on vertical axis) 
obtained from the simultaneous fit to the time evolution of RS and WS decay rates. A length scale 
is given for the axes but no axis labels and thus no central values for y'. 
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values of and y' from more data reconstructed with the most advanced tracking alignment 
constants. 

The fit function contains as one parameter the WS decay rate, defined as: 

Rws = (number of wrong-sign decays) / (number of right-sign decays) 

The simultaneous fit to the RS and WS candidate events yields the following result for the WS 
decay rate: 


Rws = ( 0.38 ± 0.04 ± 0.02) 


where the first error is statistical, the second systematic. 


This result is compatible, within its 

and 


errors, with the results obtained by E791 |n], ALEPH |^], CLEO []R|, EOCUS |H|, Belle 
is the most precise experimental measurement of Rws cnrrently available. 


4 Conclusion 


With 12.4 fb ^ of data taken by the BABAR detector in 2001, the preliminary result for y is obtained: 


y = (-1.0 ± 2.2 ± 1.7) 


(5) 


where the first error is statistical and the second systematic. This result is consistent with the 
Standard Model expectation of zero and is consistent with published values from E791 and 
EOCUS Ip and preliminary results from Belle Q and CLEO [1R|. A reduction of the statistical 
error by half is feasible in the near future by considering additional decay channels {D^ —> vr^Tr"’') 
and by including the 23 fb~^ of BABAR data collected in 2000. 

The extraction of the mixing parameters and y by measuring the time evolution of the 
wrong-sign K~^'k~) decay rate in a sample of 35 fb~^ of combined 2000 and 2001 BABAR 

data will be possible soon. 

Already with the 23 fb~^ of BABAR data collected in 2000, the most precise existing measure¬ 
ment of the wrong-sign decay rate Rws = (# decays)/RS decays) is achieved: 


Rws = (0.38 ± 0.04 (stat.) ± 0.02 (syst.)) % . 


( 6 ) 


This result is consistent with results by E791 [|^], ALEPH 
preliminary result by Belle Q. 


, CLEO m, EOCUS II and a 
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